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We establish that massive complex Abelian vector fields (mass p) can form gravitating solitons, 
when minimally coupled to Einstein’s gravity. Such Proca stars (PSs) have a stationary, everywhere 
regular and asymptotically flat geometry. The Proca field, however, possesses a harmonic time 
dependence (frequency w), realizing Wheeler’s concept of geons for an Abelian spin 1 field. We 
obtain PSs with both a spherically symmetric (static) and an axially symmetric (stationary) line 
element. The latter form a countable number of families labelled by an integer m G Z + . PSs, like 
(scalar) boson stars, carry a conserved Noether charge, and are akin to the latter in many ways. In 
particular, both types of stars exist for a limited range of frequencies and there is a maximal ADM 
mass, Mm ax , attained for an intermediate frequency. For spherically symmetric PSs (rotating PSs 
with m = 1,2,3), M maa , ~ 1.058 Mp t /p ( M max ~ 1.568, 2.337, 3.247 Mpjp), slightly larger values 
than those for (mini-)boson stars. We establish perturbative stability for a subset of solutions in 
the spherical case and anticipate a similar conclusion for fundamental modes in the rotating case. 

The discovery of PSs opens many avenues of research, reconsidering five decades of work on (scalar) 
boson stars, in particular as possible dark matter candidates. 

PACS numbers: 04.70.-s, 04.70.Bw, 03.50.-z 


Introduction. According to the latest cosmological 
data [lj, ~ 26% of the Universe’s energy content is dark 
matter (DM). The fundamental nature of DM, however, 
is unknown. A widespread viewpoint is that DM consists 
of weakly interactive massive particles (WIMPs); popu¬ 
lar candidates are the lightest supersymmetric particles, 
with masses > GeV, such as the neutralino [2j. Despite 
its success in modeling structure formation, some short¬ 
comings of the WIMPs DM model arise in small scales, 
such as the “missing satellite” and the “cuspy core” 
problems. Different proposals, some of which have been 
claimed to solve these problems (see m for reviews), in¬ 
troduce light or ultra-light bosonic particles/fields, with 
masses <C eV, which may form, gravitationally, macro¬ 
scopic Bose-Einstein condensates. Even though such 
proposals have essentially a phenomenological charac¬ 
ter, ultra-light particles (axions) are a natural ingredient 
of the Peccei-Quinn mechanism to solve the strong CP 
problem 0 and may be motivated at a fundamental level 
by string theoretical constructions - the axiverse [6j. 

Gravitationally-bound bosonic structures are thus rele¬ 
vant in the context of DM searches. The study of (scalar) 
Bose-Einstein condensates as DM candidates is often per¬ 
formed in a Newtonian limit. In the fully relativistic 
regime, such models yield gravitating solitons: (scalar) 
boson stars (SBSs). These objects, initially proposed as 
a (scalar) realisation of Wheeler’s geon idea (2,111, have 
found a variety of applications, from black hole mim- 
ickers in astrophysics, to particle models in TeV gravity 
scenarios (see [91, [k} for reviews). 

In light of recent proposals advocating massive spin 1 
particles as a DM ingredient (TTl [~Ml | , the study of the cor¬ 


responding self-gravitating structures and their dynamics 
is of special importance [151 ]. In this letter we show that 
much like massive spin 0 particles, massive spin 1 par¬ 
ticles can cluster as everywhere smooth, asymptotically 
flat lumps of energy under their own weight, producing 
gravitating solitons we dub Proca stars (PSs). Moreover, 
we observe a close parallelism between the physical prop¬ 
erties of PSs and SBSs. 

Einstein-(complex)-Proca theory. We consider one 
complex Proca field, with mass p (or equivalently, two 
real Proca fields with the same mass). It is described by 
the potential 1-form A, and field strength T = dA. We 
denote the corresponding complex conjugates by an over¬ 
bar, A and T. The minimal Einstein-(complex)-Proca 
model is described by the action: 

S = J d 4 xyf^j (yAj - \T a ^T afi - ^p 2 A a A a ^j . 


The Einstein and Proca field equations are, respectively, 
G a p = 87rGT a/ 3, \/ a P a ^ = p 2 A^, where the energy- 
momentum tensor reads: 

T a p = —F a {aFp) ~ 


+p 


A( a -A/5) 


The Proca equations imply the Lorentz condition (which 
is not a gauge choice, but a dynamical requirement): 


S7 a A a = 0 . 


The global (7(1) invariance of the action, un¬ 
der the transformation A ^ -A e la An, with a con¬ 
stant, implies the existence of a conserved 4-current, 
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j a = | [T ofi A 0 — T a ^Ap\ ■ From the Proca equation, 
V a j“ = 0. Thus a conserved Noether charge exists, Q, 
obtained integrating the temporal component of the 4- 
current on a space-like slice £: Q = f s cPxy/^gj*. 
Spherically symmetric PSs. We hist consider spher¬ 
ically symmetric solutions, taking the following line cle¬ 
ment form 

dr 2 

ds 2 = -a 2 (r)N(r)dt 2 + —— + r 2 dfl 2 , (1) 

N(r) 

where N (r) = 1—2 m(r)/r and the Proca potential ansatz 

,4 = e~ lwt [f(r)dt + ig(r)dr\ . (2) 

<r(r),m(r),f(r),g(r) are all real functions of the radial 
coordinate only and w is a real frequency parameter. 
As for SBSs the harmonic time dependence of A is cru¬ 
cial and the complex nature of A makes it compatible 
with the time-independent line element. But unlike SBSs 
there are now two independent radial functions for the 
‘matter’ held, making the analysis more involved (even 
more so in the rotating case). The Proca equations yield 

d_ ( r 2 [f - wg] 
dr 1 a 


h 2 r 2 f 

aN 


wg- f = 


y 2 o 2 Ng 


where ' denotes radial derivative. These two equations 
imply the Lorentz condition constraint (which determines 
f(r) in terms of the remaining functions). The essential 
Einstein equations read 

1 


m = 47rGr 2 


— = AnGry 2 
a 


if - wg) : ± 2 
2 a 2 2 M 


g 2 N- 


f 2 

Na 2 


9 


CA 


N 2 a 2 J ' 

For the ansatz 0-©, the Noether charge reads 
Q = / 0 °° dr r 2 g(r) 2 a(r)N(r), and the energy den¬ 

sity measured by a static observer is 


_ T ‘ - if ~ W9) 2 + -a 2 
“ 2a 2 + 2 M 


g 2 N- 


P\ 


Na 2 J 


Finally, PSs satisfy the virial relation: 


dr r 2 a 


M 9 


/ 2 (41V — 1) \ (wg-f ) 2 


t 2 N 2 


= 0, 


used to test the numerical accuracy of the results. This 
relation can also be used to rule out non-gravitating so¬ 
lutions, i.e ‘Proca-balls’ without backreaction 0. 
Asymptotic behaviours. An analysis of the held equa¬ 
tions both near the origin and at spatial inhnity reveals 
smooth behaviours. Close to r = 0 we hnd 


fo 


w 


f(r) = fo + ^^ 2 -^)r 2 +0(r*), 

9( r ) = ~i^rr + 0(r 3 ) , 

6(7 0 


, ^ 4ttG/ 0 V 3 
m(r) = - —r 


6 ct 2 


G(r 5 ) 


c(r) = cr 0 + 


2 . .2 


47T Gf 2 y 
2cto 


■r 2 + G(r 4 ) , 


where fo, cr 0 are constants, the values of f(r) and a (r) 
at the origin, respectively. As r —> oo, we hnd 


f(r) = c 0 - 


g(r) = c 0 


o-ry/ll 2 -W 2 


„-ry/ ii 2 -w 2 


\Jp 2 — W 2 r 

m{r) = M + ..., 

2 2 -2 ry/^-w 2 

log^) = -4^ 2( „G^ ) 3 / , -7-+ •■ 


where M is the ADM mass and Co is a constant; observe 
that w < y, which is a bound state condition. 
Numerical results. The solutions that smoothly inter¬ 
polate between the two above asymptotic behaviours are 
found numerically. In numerics - and in the following - 
we set y = 1, 47tG = 1, by using a scaled radial coordi¬ 
nate r —> rp, (together with m —> my, w —> w/y ) and 
scaled potentials / —> f\/AirG, g —> gy/AnG. The equa¬ 
tions are solved by using a standard Runge-Kutta ODE 
solver and implementing a shooting method in terms of 
the parameter /(0). 

In Fig. |T]we plot the ADM mass, M, and the Noether 
charge, Q, as a function of the scalar held frequency, w. 
As w —> 1, the mass and Noether charge of the solutions 



FIG. 1. ADM mass M and the Noether charge Q for spher¬ 
ical PSs vs. the Proca held frequency w. The inset shows M 
vs. the shooting parameter /(0), which is analogous to the 
curve for M vs. the central value of the scalar held for SBSs. 

vanish, but Q/M —> 1. Approaching this limit, PSs - 
like SBSs - become spatially diluted with an effective 
size much larger than their Schwarzschild radii 0 , and 
trivialize at w = 1. Reducing w from this maximal value, 
PSs become more compact; both M and Q follow a spiral, 
towards a central critical configuration located around 
w — 0.891. The maximum of both M and Q occurs for 
Wmax — 0.875, with M max ~ 1.058 < Q m ax — 1.088. 
We observe that M < Q from w = 1 down to almost the 
minimal allowed frequency. In the lower part of the spiral 
M > Q, and the binding energy 1 — M/Q is negative, 
corresponding to a region of “energy excess”, where the 
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solutions must be unstable against perturbations. All 
these features follow closely those of spherical SBSs. 

We remark that for the family of solutions exhibited 
in Fig. ffl /( r) has one node JT^] and g(r ) is nodeless. 
These are the fundamental solutions. There are also ex¬ 
cited solutions, with more nodes for these functions, but 
we did not attempt to study them in any systematic way. 
Studies of SBSs suggest that excited solutions are unsta¬ 
ble [19]. In Fig. E| we show the behaviour of physical 
and profile functions for two PS solutions. In particular 
observe the energy concentration near r = 0. 


allows to write the perturbed metric as: 

ds 2 =—cr 2 (r)N(r) [l - eh 0 (r)e~ int ] dt 2 

r 2 dfl 2 


dr 2 [l + ehi{r)e~ iM ] , 2 ,_ 2 . 


N(r) 

the vector field is perturbed as 

-ifit e/iO) -Me/ 2 (r) 


A =e 


f(r) 


dt 


ig{r) 


j^im e 9i(r) + ieg 2 (r) 


dr 


(3) 


(4) 




FIG. 2. Energy and Noether charge densities, electric field 
(E r = | Frt |) and profile functions (inset) of two PSs solutions 
with M = 0.745, w = 0.817, /(0) = 0.394 (top panel) and 
M = 1.016, w = 0.839, /(0) = 0.165 (lower panel). 


Stability. A positive binding energy - as that observed 
along the upper branch of Fig. |T] does not guarantee 
linear stability. For SBSs, a perturbative stability anal¬ 
ysis shows that the solutions are only stable from the 
maximal frequency w = 1 down until the frequency at 
which the maximal mass is attained, Wm,w Beyond this 
point an unstable mode develops |20| (2l| • An analogous 
conclusion holds for PSs, as we now show by considering 
their linear radial perturbations. 

We assume that all perturbations have a harmonic time 
dependence of the form e~ lQt , with 12 being the charac¬ 
teristic vibrational frequencies of the PS. Gauge freedom 


where ho, hi, fi, / 2 , g\ and g 2 are radial perturbations 
around the background solution, and e is a small, book¬ 
keeping parameter. The perturbed Einstein-Proca sys¬ 
tem can be reduced to a pair of second order ODE’s. Im¬ 
posing regularity of the perturbations at the origin and at 
infinity, the resulting system is a two dimensional eigen¬ 
value problem for 17 and one other constant which we 
have chosen to be the value of ho at the origin. A numer¬ 
ical solution is then obtained by a two dimensional shoot¬ 
ing, with the result shown in Fig. [Sj where we plot I2 2 ver¬ 
sus the PS’s ADM mass around the maximal mass. For 



FIG. 3. Fourier frequency 2 of the perturbations as a func¬ 
tions of the PS’s total mass. The critical mass at which PSs 
become unstable corresponds to the star’s maximum mass. 


f2 2 > 0 the mode is purely real and corresponds to stable 
normal modes. This occurs for w > w max . For I2 2 < 0 
the mode is a purely positive imaginary number, indicat¬ 
ing that these configurations are unstable (c/. Eqs. @ 
and (^). Thus, in complete agreement with the SBSs 
case 20|-|23|. we can see that the star’s maximum mass 
corresponds to a transition point between stable and un¬ 
stable configurations. For unstable configurations with 
mass 0.1% below the threshold, the instability timescale 
r is already smaller than r < 100M. Thus, we expect the 
unstable branch of PS to share similar properties to those 
observed in the scalar case: configurations which reach 
this branch will either quickly collapse to black holes or 
migrate back to the stable branch via mass ejection (the 
“gravitational cooling” mechanism) [22], IHHUl • 
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Rotating PSs. We now turn to rotating PSs. These 
solutions are found with an axi-symmetric metric ansatz 

ds 2 = -e 2 F °^dt 2 + e 2 F ^ e \dr 2 + r 2 d9 2 ) 

+e 2 F ^r 2 sin 2 d (dip - ^^-dtj ' . 

The Proca held ansatz is given in terms of another four 
functions ( Hi,V) which depend also on r, 9 

A = (^ydr + H 2 d9 + iH 3 sin 9dp + iVdt)j e dm V -wt), 

with m G Z + . The Einstein-Proca equations are solved 
with the following boundary conditions, which we have 
found to be compatible with an approximate construc¬ 
tion of the solutions on the boundary of the domain of 
integration: d r Fi\ r _ Q = W | r=Q = i?i| r =o = U| r=0 = 0; 
also = Hi\ r —oo U| r=00 = 0; and 

deFi\ e=Q = dgW | 0=o = ffi\e=o,n = V\g= 0 ^ = 0 (note 
that for m = 1 the conditions for H 2 ,H 3 are different, 
with dgH 2 \ g _ Q ^ = dgH 3 \ g _ Q ^ = 0). All rotating PSs 
we have constructed so far are symmetric w.r.t. a reflec¬ 
tion along the equatorial plane. Odd-parity composite 
configurations, however, are also likely to exist. 

As usual, the ADM mass M and angular momentum 
J are read off from the asymptotic expansion, g tt = —1 + 

2 GM/r + ..., g v t = —2GJ sin 2 9/r + _ In analogy 

with the SBSs case (13], one can show that the angular 
momentum is a multiple of the Noether charge, J = mQ, 
but in contrast to the SBSs case, the angular momentum 
density T * and the Noether charge density j* are not 
proportional any longer. 

The Einstein-Proca equations for this rotating case are 
quite involved and shall not be presented here. They are 
solved numerically, subject to the above boundary con¬ 
ditions, by using the elliptic PDE solver FIDISOL [28j 
based on the Newton-Raphson procedure, employing a 
compactified radial coordinate x = r/(l + r) [ 23 . In 
Fig. U we exhibit the mass vs. angular momentum, to¬ 
gether with the mass and the Noether charge vs. the 
Proca field frequency, for rotating PSs with m = 1,2,3 
(inset). Both panels exhibit the familiar shape from the 
study of SBSs. In particular M and J are always pos¬ 
itively correlated [3(|. The physical quantities for the 
largest mass PSs are exhibited in the following table, for 
m = 0 (spherical) and m = 1, 2,3 (rotating), where they 
are compared with the results in the literature for the 
corresponding largest mass SBSs (last column); one ob¬ 
serves that M max is always smaller for the latter. 



U)max 

M-max 

Qmax 

Jmax 

Mmax SBSs 

m = 0 

0.875 

1.058 

1.088 

0 

0.633 [2 

m = 1 

0.839 

1.568 

1.626 

1.626 

1.315 [27, 31] 

m = 2 

0.767 

2.337 

2.461 

4.921 

2.216 [31] 

m = 3 

0.667 

3.247 

3.489 

10.467 

not reported 



FIG. 4. M vs. J diagram for rotating PSs with m = 1,2,3. 
The inset shows M (solid red lines) and Q (dashed green lines) 
vs. w, where m = 0 corresponds to the spherical PSs in Fig. [T] 

Other physical properties of the rotating PSs mimic 
those of SBS: the rotating PSs solutions spatially delocal¬ 
ize as w —> 1 (similarly to the static case) and trivialize in 
that limit; the energy density has an essentially toroidal 
distribution; also, for m = 1, neither the energy density 
nor the Noether charge density vanish on the symmetry 
axis [ 32 I ] . 

We have not attempted to study in detail the stability 
of rotating PSs. For rotating SBSs, catastrophe theory 
arguments [33| support a similar conclusion to that ob¬ 
tained in the static case: the solutions are stable from 
the maximal frequency down to the point where the max¬ 
imal mass is attained. We expect a similar conclusion to 
hold for rotating PSs. In the rotating case there can also 
occur ergoregion instabilities Q; for SBSs ergoregions 
only appear in the unstable branch of solutions and an 
analogous situation is true for PSs (35| . 

Discussion and Outlook. The discovery of PSs sug¬ 
gests various novel directions of research, besides the de¬ 
tailed analysis of the solutions presented herein; let us 
mention two. 1) In the scalar case, a single (rather than 
complex) real field allows the existence of very long lived 
quasi-solitons - oscillatons ; the same is true for the Proca 
case [26]. 2) Can one add a black hole inside a PS, like 
for other gravitating solitons? While for spherically sym¬ 
metric configurations one can show the absence of black 
hole solutions, this is possible inside a spinning PSs, again 
in complete analogy with the SBSs case [36l [381 . 

A final remark relating PSs with DM phenomenology. 
Since ordinary matter, which constitutes only ~ 5% of 
the Universe’s energy content, is comprised of over ten 
elementary particles, it is reasonable to admit that DM is 
composed of different kinds of fundamental entities, but 
which, when gravitationally clustered into macroscopic 
lumps, display some universality. If Bose-Einstein con¬ 
densates of ultra-light scalar fields are viable DM models, 
the existence of PSs akin to SBSs suggests the former may 
be another DM component. In this context it would be 
interesting to study the Newtonian limit of PSs. 
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